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Quantum fidelity between two density matrices, F(p\,p2) is usually denned as the trace of the 
operator T — ^/- s /pip2 y/pi- We study the logarithmic spectrum of this operator, which we denote 
by fidelity spectrum, in the cases of the XX spin chain in a magnetic field, a magnetic impurity 
inserted in a conventional superconductor and a bulk superconductor at finite temperature. When 
the density matrices are equal, pi = pi, the fidelity spectrum reduces to the entanglement spectrum. 
We find that the fidelity spectrum can be a useful tool in giving a detailed characterization of different 
phases of many-body quantum systems. 



I. INTRODUCTION 



A quantum system in a pure state is described by a density matrix which is just a projector onto that state. At 
zero temperature it is the projector to the groundstate of the system. In general, the Hamiltonian of the system is a 
function of some parameters which determine the groundstate. The quantum fidelity between two states (for two sets 
of parameters) is, in this simple case, the absolute value of the overlap between the groundstates for the two sets of 
parameters. When the system is in a mixed state the density matrix is more complex. Typical situations that lead to 
mixed states are: i) reduced density matrices where a trace over some degrees of freedom is carried out, or ii) systems 
at finite temperatures where the density matrix may be taken as the Boltzmann factor over the energy eigenstates. 

In general the quantum fidelity [lj between two states characterized by two density matrices p\ and p2 may be 
defined as the trace of the fidelity operator, T , 



C / 

C3 F(pi,p 2 ) = TrJ 7 = Tryjy/ftp2y/pi- (1) 

One may also consider the spectrum of the fidelity operator F(pi,p2)- Its set of eigenvalues A^, which we denote 

fidelity operator spectrum, and — In A^, which we call the fidelity spectrum, may provide more information as compared 

^H to the fidelity (its trace), in a way parallel to the extra information provided by the entanglement spectrum [2], as 

compared to the von Neumann entropy. 

In the case of pure states p\ = \GSi)(GSi\ and p2 = \GS2)(GS2\ the fidelity is just the norm of the overlap 
CO 

ON F(p 1 ,p 2 ) = |(G5i|G5 2 >|. (2) 

[^. If the two states are the same it is just the normalization of the state (taken as 1). In the case of two equal mixed 
^^ states the fidelity is just 

^H 

^H F(p,p) = Tip=l (3) 



and the operator T in this case has a set of eigenvalues, A^ = Af, such that — In A^ is called the entanglement spectrum, 
and has received considerable attention lately |2J. 

In this work we will analyse the fidelity spectrum for several physical systems paying particular attention to the 
vicinity of quantum phase transitions (QPT) as well as the properties characterizing their quantum phases. 

The interplay between quantum information and condensed matter physics has been extensively considered using 
entanglement as a measure for the behavior of many body systems |3J. The distinguishability between states has 
been used as a possible criterion to study quantum phase transitions |4J. By its own nature, fidelity between pure 
groundstates signals a change of state as one approaches a quantum phase transition |5J . The fidelity between mixed 
states has also been used as a signature of quantum phase transitions [6, 7 J and to distinguish between different states 
of matter at finite temperatures [8 J . A standard measure of entanglement in a system is the von Neumann entropy. 
However, as argued in Ref. |2J, more information about a mixed state is obtained if the entanglement spectrum is 
analysed. Considering reduced density matrices where part of the degrees of freedom are integrated over, such as 
dividing the system in real space into two parts A and B, it was shown in the context of the quantum Hall effect 
[2j [9] and in the context of coupled spin chains [10J, that the groundstate entanglement spectrum of A contains 
information about excited energy states of the frontier of the subsystem A. In particular, in the quantum Hall effect 
the entanglement spectrum of the bulk system has a low- lying structure of levels that matches the edge states, and in 
the case of two coupled Heisenberg spin chains, considering the subsystem A as one of the chains, the entanglement 



spectrum has a structure that matches the energy excitations of a single Heisenberg chain. Other partitionings of the 
system have been proposed that lead to further information [11 J and considering a partitioning in momentum space 
it was shown that information about energy excitations of a single Heisenberg chain is contained in the groundstate 
wave function through the entanglement spectrum [12J. 

In this work we consider the fidelity spectrum of various systems. While the entanglement spectrum has some 
relation to the energy spectrum of the edge states or even bulk states, the fidelity spectrum contains information 
about which eigenvalues have a larger contribution to the distinguishability between quantum states. We start by 
considering two systems at zero temperature: a magnetic impurity inserted in a conventional superconductor and the 
XX chain in a magnetic field. In the first case, described in Section II, the magnetic impurity is coupled through 
a spin interaction to the spin density of the conduction electrons, tuned by a coupling J. As previously discussed 
[7J[T3j[14], as the coupling J grows, the system goes through a first order phase transition. At this point the system 
becomes magnetized and various quantities such as local density of states, spin content, gap function and quantum 
information measures can be used to detect this transition. For instance, various entanglement measures [15J and the 
partial state fidelity itself |7J have been used before. In the case of the XX chain describing spins 1/2 confined to 
a plane (xy) and with a transverse magnetic field h, aligned along the z direction (Section III), there is a quantum 
phase transition from a XX phase, where the spins are aligned in the xy plane, if the magnetic field is small, and an 
Ising-like phase where the spins point along the field direction, if the field is strong enough. Considering the coupling 
between the spins as the energy scale, the second order transition occurs at the point h c = 1. This transition is also 
signaled in various ways such as the decrease of the fidelity near the critical point |4J. In Section IV we consider 
thermal states in a conventional superconductor. 

II. MAGNETIC IMPURITY IN A SUPERCONDUCTOR 

Consider first a classical spin immersed in a two-dimensional s-wave conventional superconductor. We use a lattice 
description of the system. In the center of the system, i — l c — (0,0), we place a classical spin along the z-direction 
S = Se z , with no loss of generality. The Hamiltonian of the system is given by 

H = -J2 *tf c L c i* - £ f J2 C *W + Yl ( A *4r c l + A i c a c n) - J2 Jc la aZ ^ c ic^ , (4) 

(ij)a i<J i <j<j' 

where the first term describes the hopping of electrons between different sites on the lattice, Sp is the chemical 
potential, the third term is the superconducting s-pairing with the site-dependent order parameter A^, and the last 
term, with J > 0, is the exchange interaction between an electron at site i = l c and the magnetic impurity. The 
hopping matrix is given by tij = tSj^+s where S is a vector to a nearest-neighbor site. Note that both indices 
z, j G {1,2,... N} specify sites on a two-dimensional system (N is the number of sites). We take energy units in terms 
of t (t = 1), and choose £p = — 1. 

If we divide the whole system in two subsystems, say A and £?, then the partial mixed state, given by the reduced 
density operator pA for the subsystem A, is defined as 

Pa = Tr B p, (5) 

where Tr#[-] represents the partial trace evaluated over the Hilbert space %b of the subsystem B. We take A to be 
one site, either the impurity site or an arbitrary site in the bulk far from the impurity such that p\ and p2 of Eq. ([1]) 
are one-site density matrices. 

In many-body systems second quantization is the natural way to perform any calculation. The matrix elements of 
the density matrix are simply defined in terms of correlation functions of the whole system. For instance, in the case 
of the single-site partial mixed states using local basis states B = {|0), | tt)? It)? It)}? which denote the four possible 
states — unoccupied, double occupied, single occupied with an electron with spin up and single occupied with an 
electron with spin down, respectively — it can be shown that the corresponding density matrix reads as [7] 
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where the index i denotes the site. The spin and the charge parts decouple. The spin part couples the two spin 
orientations (single occupied states) and the charge part couples the empty and doubly occupied states. The diagonal 
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Figure 1: (color online) Fidelity operator spectrum at the impurity site (left) and at a bulk site (right) as a function of the 
spin coupling, J, where one density matrix is calculated at J and the other at J + 5 J, where 5 J = 0.05. In black and red is 
the charge contribution and in green and blue the spin contribution. 
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Figure 2: (color online) Fidelity operator spectrum: charge eigenvalues in black and red; spin eigenvalues in green and blue. 
System size is 15 x 15. a) One of the sites is at the impurity, i± = l c , and the other site, %2 : is arbitrary. In the left panel 
J\ — Ji — 1.5, and in the right panel J\ — J2 = 2.5. b) One of the sites, ii, is a site in the bulk and the other site, 22, is 
arbitrary. The other parameters are the same as in a). Recall that J — 1.5 < J c and J — 2.5 > J c . 



terms of the matrix describe the number of empty sites, the number of doubly occupied sites, the number of spin up 
sites and the number of spin down sites, respectively. The sum of the diagonal terms is equal to 1 due to normalization. 
The matrix is easily diagonalized and the fidelity between two different one-site states obtained straightforwardly. 

We will consider two density matrices of the form pi{Ji]i\), and p2{J2'^2) for two values, J\ and J2, of the spin 
coupling J, and for two sites, i\ and %2- 

In Fig. [T] we show the fidelity operator spectrum as a function of the coupling between the magnetic impurity and 
the electronic spin density taking J\ = J, J<i = J + 0.05 where i\ = 12 = i is the impurity site or a bulk site. The 
charge and spin parts separate and in the charge part there is the empty and the doubly-occupied contributions and 
in the spin part the spin up and spin down contributions. As the coupling J grows there is a discontinuity in the 
eigenvalues. This discontinuity is associated with the quantum phase transition previously discussed. The sum of the 




Figure 3: (color online) Fidelity operator spectrum as a function of space. System size is 15 x 15. Top panels J — 1.5. Lower 
panels J = 2.5. From left to right, the first two panels are the charge eigenvalues and the right two panels the spin eigenvalues, 
for lattice sites i\ — l c and 12 an arbitrary site. 

four eigenvalues is the fidelity, as discussed before. There is a discontinuity at the QPT both in the total fidelity, in 
the charge and spin parts and in the individual eigenvalues as well. As we can see, the discontinuities occur mainly 
in one of the charge eigenvalues and in one of the spin eigenvalues. As one crosses the QPT to a regime where the 
impurity captures one electron breaking a Cooper pair we see that the main contribution to the discontinuity in the 
charge part comes from the doubly-occupied states, where there is a significant increase at the QPT. In the same way 
there is a significant decrease in the spin up eigenvalue, leading to a smaller spin contribution beyond the QPT. For 
small values of J there is a screening of the perturbation induced by the magnetic impurity in the superconductor. A 
small fidelity means a higher degree of distinguishability. For small values of J the spin down eigenvalue is small but 
the spin up is still high. Beyond the QPT both contributions are small. This indicates that the transition is mainly 
of spin character: as the fidelity tends asymptotically to 1 away from the QPT and the spin eigenvalues are small, 
the charge eigenvalues have to compensate, however, mainly through the doubly-occupied contribution, as expected. 
Note that, in the bulk, the transition is quite small. The physics is very local, centered around the impurity site. 

In Fig. [2] we consider that the coupling is fixed, but the two density matrices are calculated at different sites i\ and 
^2' pi(J,h) and ^2(^,^2)- The various cases are specified in the caption of the figure. In these figures the horizontal 
axis is the lattice site. For each site there are four eigenvalues. The central point is the impurity site. The symmetrical 
peaks close to the central point are the neighbors of the central point (please note that here we number the lattice 
sites of the 2d system as a Id system sequentially row by row; so the lattice nearest neighbors in the x-direction are 
close by neighbors but the neighbors in the y direction are far apart). The fluctuations around the central peak are 
of course better seen in a 3d plot. This is shown in Fig. [3] where we compare the four eigenvalues with the results of 
Fig. [TJl. As discussed above the empty site and the spin down contributions are small and do not change much as we 
change from J = 1.5 (below the QPT) to J = 2.5 (above the QPT). 



III. XX SPIN-1/2 CHAIN IN A TRANSVERSE FIELD 



The XX spin 1/2 model has been solved exactly via the Jordan- Wigner transformation where it is reduced to a 
system of free spinless fermions [16J. The correlation functions were also calculated [17J as well as the reduced density 
matrix of a system of L contiguous spins [18-20J. The information theoretic approach in terms fidelity, Fisher metric 
and Chernoff bound was applied to the XY model in [^ [^TTBfj] . 

The Hamiltonian we will consider here is of the form 
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(7) 



Here the spin operators are described by Pauli matrices, h is the transverse magnetic field and 7 is the anisotropy. 
We will simplify and consider 7 = and we take < h < 1. We will consider a block of L contiguous spins. The 



reduced density matrix of the block can be written as [18J 



p A = n (^bh + ^b. 



(8) 



The operators bi are spinless fermionic operators with a 2-state space (the eigenvalues of the number operator U{ = b\bi 
are 1 and 0, corresponding to occupied or empty state, respectively). The eigenvalues of the reduced density matrix 
are 2 L in number as the result of the direct product of the i = 1, • • • , L subspaces. Defining two Majorana operators 
of the form 
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in terms of the spin operators, it has been shown that 

(GS\c m c n \GS) = 5 m , n + i(B L ) r 



The matrix Bl is written as 



with 
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We have that g\ = g_^, go = (2ip c /ir) — 1 and gi^o = (2 /lir) sin ltp c , where cp c = arccos(/i). Defining new Majorana 
fermions through the transformation d = Vc and imposing that 



with 
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B L = VBjV T = Q < 



\ / ? 



o 1 
-1 o 



(13) 



(14) 



where Q is a diagonal matrix with L diagonal elements v\, leads to the diagonal form Eq. (|sj) of the reduced density 
matrix, p^, having defined L complex fermionic fields like 



bi 



id 



2Z+1 



(15) 



The transformation V that block diagonalizes the problem depends on the Hamiltonian parameters through the 
numbers g\. 

In order to calculate the fidelity operator T = ^/ sfpi p2 \fpi one needs to consider the product of two reduced density 
matrices for two different magnetic fields. Even though each of the reduced density matrix can be diagonalized, the 
transformations needed to diagonalize p\ and p2 are different. The diagonalization of p\ (p2) is obtained introducing 
a matrix V\ (V2). To obtain the spectrum of T we rewrite the expression for the diagonalized reduced density matrix 
p2 in terms of the fermionic operators of the matrix p\ which leads to p2 = e~ H<2 , where 
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(16) 
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Figure 4: (color online) Fidelity spectrum — lnA^, for system size L — 6. Left: entanglement spectrum for different values 
of magnetic field. From bottom to top h = 0.6, 0.7, 0.8, 0.9, 0.95, 0.99. Middle: fidelity spectrum for Sh = 0.01 and for 
h = 0.7, 0.8, 0.9, 0.95, 0.99. Right: fidelity spectrum for pairs of the magnetic field as hi = 0.5, h 2 = 0.5, 0.6, 0.7, 0.8, 0.9. 
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where we defined V by 
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(18) 



The diagonal elements 1/2,1 are the diagonal elements of Cl for the magnetic field of p2- Diagonalizing i^2, we obtain 
p2, expressed in the eigenbasis of pi, and obtain the spectrum of J 7 , as intended. Note that this requires diagonalizing 
a 2 L x 2 L matrix which is much larger than the L x L matrix required for the entanglement entropy. 

In Fig. [4] we show the fidelity spectrum — In A^. In the left panel we consider p\ = p2 and the operator T is just the 
density matrix. Therefore, its logarithmic spectrum is the entanglement spectrum. In the middle panel we consider 
the fidelity spectrum and in the right panel we compare the entanglement spectrum with the fidelity spectrum. As 
discussed before [12J there is no clear structure in the entanglement spectrum since we are considering a real space 
block of spins. The fidelity spectrum when we consider two very close values of the magnetic field is very similar to the 
entanglement spectrum. In both cases as we increase the magnetic field and approach the critical point that separates 
the XX phase from the Ising phase, the entanglement or fidelity spectrum increases considerably. This implies that 
the fidelity operator spectrum is decreasing fast. The same can be seen when we compare the fidelity spectrum with 
the entanglement spectrum in the right panel. Here we are considering two values of the magnetic field h\ and hi 
associated respectively with p\ and p2 that differ by a finite amount. Therefore there is a significant difference as the 
difference between h\ and Y12 increases, even though we are far from the critical regime. The further the two points 
are the smaller the fidelity and the fidelity operator eigenvalues, A^, should be. 

In order to analyse the spectrum we calculate the moments of the distribution of the eigenvalues [27J defined as 



M n = J2K- 



(19) 



The moment of order n = 1 is the fidelity. We also calculate the von Neumann entropy of the fidelity 



and the Renyi entropies 
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Note that in the case of the entanglement (pi = P2) the von Neumann entropy can also be obtained as 
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Figure 5: (color online) First five moments of the entanglement operator (density matrix) and fidelity operator spectra for 
blocks of sizes L — 1, 2, 4, 6, as a function of magnetic field. In the top left panel we consider the moments of the entanglement 
operator spectrum. In the top right panel we consider the moments of the fidelity operator spectrum for two close- by values of 
h. In the bottom left (right) panel we consider the fidelity operator spectra for hi = 0.5 (hi = 0.98) and hi arbitrary. 



which only involves a sum with L terms. 

In Fig. [5] we present the first five moments for several choices of pairs of magnetic fields (hi ^2) and for several 
block sizes, L. In the case of the entanglement spectrum, when hi = /i 2 , the first moment is just the trace of the 
density matrix (top left panel). When hi 7^ fo 2 it is the fidelity and therefore it is very close to 1 except in the vicinity 
of the phase transition where we consider /i 2 — hi — 5h,5h = 0.01 (top right panel). The higher momenta are very 
similar between the entanglement spectrum and the fidelity spectrum, in this case. A detailed study of the moments 
and Renyi entropies has been carried out in Ref. [28J. Considering two values of the magnetic fields further away 
from each other the moments change in structure. Fixing for instance hi = 0.5 (bottom left panel) or hi = 0.98 
(bottom right panel) and varying ft, 2 we find that when ft 2 crosses hi there is a sharp increase which equals one for the 
first moment (trivially since for this case the first moment is the trace of the density matrix). This maximum is also 
observed in the higher moments of the spectrum. A similar information can be obtained from the Renyi entropies. 
As explained above, Si is just the von Neumann entropy and the other entropies are proportional to the logarithms 
of the moments of the spectrum. In Fig. [6] we compare various cases for L = 1,6. As the magnetic field approaches 
the critical point we find that the various entropies have a minimum near the critical point showing that the partial 
state density matrix signals the QPT, as previously obtained. In the lower panels we consider the Renyi entropies 
associated with the fidelity operator. The structure is in general more complicated. Fixing hi = 0.98 close to the 
critical point, as hi approaches hi the entropies tend to those corresponding to the entanglement spectrum. Far from 
this point the entropies differ considerably and depend strongly on the block size. For instance for L = 1 they are of 
the order of the entropies for the entanglement spectrum but for L = 6 the entropies are considerably higher, except 
the von Neumann entropy. Also, note that the Renyi entropy Si in the fidelity case has a depression at small magnetic 
fields. 
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Figure 6: (color online) Renyi and von Neumann entropies for L — 1 (left) and L — 6 (right), for the first five moments, as 
a function of magnetic field. In the top panels we consider the entanglement operator spectrum and in the bottom panels the 
fidelity operator spectrum as a function of /12, for h\ — 0.98. Note the decrease of the entropies close to the quantum critical 
point, both for the entanglement case and the fidelity case. 
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Figure 7: (color online) Block fidelity susceptibility as a function of magnetic field for L = 1, 
increase near the critical point. 



. . , 6 which shows the sharp 



We also calculate the fidelity susceptibility, introduced in [29J and [23J, with its geometrical meaning discussed in 
[231 1301 • Global fidelity susceptibility for the XY model was discussed in [UH3]. In the present paper was calculate 
the fidelity susceptibility associated to reduced density matrices of blocks of spins. It is defined as 
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In Fig. [7] we present the fidelity susceptibility associated with blocks of different sizes as a function of the magnetic 
field. The divergence of the susceptibility is clearly seen as we approach the QPT. 
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Figure 8: (color online) Exponential entanglement spectrum for a bulk superconductor as a function of momentum, labeled 
sequentially row by row in the Brillouin zone, in the normal phase (A = 0, left panel), and in the superconducting phase 
(A/0, right panel), at temperature T. Note that in the right panel the lowest charge eigenvalue is smaller than the spin 
eigenvalues. 
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Figure 9: (color online) Fidelity operator spectrum for a bulk superconductor as a function of momenta, between the normal 
phase and the superconducting phase, at the same temperature (left) and at different temperatures (right). 



IV. THERMAL STATES OF A CONVENTIONAL BCS SUPERCONDUCTOR 



In this section we consider a conventional s-wave superconductor at finite temperature described by the effective 
mean-field BCS Hamiltonian 
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(24) 



with Afc = — V(c-ki c kt)i wnere the lattice-mediated pairing interaction is constant and non- vanishing between elec- 
trons around the Fermi level only. The density matrix is given by |8J 
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where T is the temperature, h k = (ft+ ftj^/ig) = (2A%/T,2A k /T,-2e h /T), f k = (T+,T-,T fc °), K 
A^fr/e) and e k = £ k —p- The norms of the vectors h k are given by h k = 2E k /T : with E k 
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hk = hk(T,V) are functions of both the coupling constant V and the temperature T, through the gap parameters 
Afc = Afe(T, V) and the chemical potential fi. By n ka = c\ a c k o- we denote the one-particle number operators, while 
by b\ = 4t c -/c; and b k = C- k iC k ^. The T k operators are given by T^ = b\, T fe " = b k and 2T^ + 1 = (n k ^ + n_^) and 
form a su(2) algebra. 
The fidelity is given by 



Twrn p^ L TkJ k T kp ^f k \i/2] 

F(p a ,p b ) = Tr[(py 2 ^ /2 ) 1/2 ] = [(iU J 1\ } l = Tr 

KH J LV a a ) J n fe (Tr[e^^]Tr[e^^])V2 



ri(^) 



1/2 



(26) 



10 




Figure 10: (color online) Top (bottom) panels: normal (superconducting) phase exponential entanglement spectrum. In the 
three panels we show three eigenvalues since the two spin eigenvalues are degenerate. 




Figure 11: (color online) Fidelity operator spectrum between superconducting and normal phases. In the first 3 panels we 
show the 3 eigenvalues and in the last one the fidelity (smaller than 1 close to the Fermi surface). Note that the temperature 
is the same, with A either finite (superconducting phase) or zero (normal phase). 



As for every k the operators X^ form a su(2) algebra, and therefore by exponentiation define a Lie group, we can write 
e ^rT k e b k T k e ^T k _ e 2c fc T fc p or eac j 1 va l ue f the momentum we get a 4-dimensional space respecting to momentum 
states that are empty, doubly occupied or singly occupied by a spin up or a spin down electron. The space is therefore 
of the type B = {|0), | tl), 1 1)> I i)}> similarly to the problem of the impurity in a superconductor. The fidelity operator 
is then easily diagonalized in this 4x4 subspace. We will study the possible eigenvalues of the fidelity operator for 
each momentum. As before, the charge and spin parts separate. Moreover, in this problem the two spin components 
are degenerate. Therefore it is enough to look at three eigenvalues (two for the charge part and one for the spin part). 
The results are presented in the Figs. 8-12 and discussed in more detail below. _ 

The fidelity operator for each momentum value (denoted k-fidelity operator) is of the form T^ = VAk, where 





/ a-k 


fa o 


1 


Pk 


7fe 


3fe 


, o 


1 




V o 


1 



(27) 
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The matrix elements are given by 



&k 



cosh El jT cosh E\ jT + sinh E% jT sinh E\ jT 



A£A 



k^k 



etc 



k^k 



E a k E b k 



+ sinh E a k /T cosh E b /T^ + sinh E b /T^ 



E a k 
(cosh E^/T - l)smhE b k /T % k 



pa pb pa 
t k t k t k 



E a k E b k 



El 



7fe 



cosh El/T cosh E b k /T + sinh E%/T sinh ££/T 



- sinh El/T cosh E b k /T^ - sinh E\/T^\ 



A£A£. 



6 fc 6 fc 



^ 



(cosh££/T - l)sinh^/T A ^ Afe 



El 

pa pb pa 
t k t k t k 



E a k E b k 



El 



fa 



sinh E a k /T cosh E b k /T^ + sinh E b /T^ 



At 



+ (cosh^/T-l)sinh^/T 



A£A£ + elel A a 



k 

fc c fe A fe 



^k 



^ a 



and 



L> fc = 2 (1 + cosh El/T) 2(1 + cosh ^/T) . 

The eigenvalues of the fc-fidelity operator, T k are therefore of the form 

(j)% 0\ 

1 rt 

\fD~ k 10' 

\ 1/ 



where 



V± 



(a k + 7fe ) ± ^(afc - 7fe) 2 + 4/9, 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



We will be interested in situations where pi and p2 correspond to points in parameter space, which we choose to be 
the temperature, T, and the gap function, A&, that are far apart and may be in the same or different thermodynamic 
phases. 

In Fig. [8] we present the /c-fidelity operator spectrum (of the operator F^) for the case when the two density matrices 
are equal (^'exponential entanglement spectrum") where we compare the system in the normal phase (left) with the 
superconducting phase (right). The horizontal axis is an index over the eigenvalues and in the vertical axis we plot 
Afc. We only plot three eigenvalues because the spin eigenvalues are degenerate. For each label the sum over the four 
eigenvalues is 1 due to normalization. In both phases the higher eigenvalue is the charge eigenvalue corresponding 
to empty sites (this will be discussed later on). In the normal phase the lowest eigenvalues merge into the higher 
eigenvalues but in the superconducting phase the energy gap is clearly visible. Note that the eigenvalues are now 
labeled by the momentum. There is no partitioning of the system in real space but there is a partitioning of the 
system in momentum space (since the system can be block diagonalized) . Recall however that here the mixed state 
originates in the thermal states. We stress that we are not plotting the fidelity operator, T = EL ^~fc> eigenvalues. 
In Fig. [9] we consider two different density matrices where we plot the fc-fidelity operator spectrum where one of the 
density matrices corresponds to a point in phase space in the normal phase and the other in the superconducting 
phase. In the left panel the temperatures, Ta^T^, are the same (this can be obtained for instance considering two 
coupling constants) and in the right panel the temperatures, Ta^T^, are different. In the first case the gap is still 
clearly visible. When the temperatures are different the gap remains the same. There is a small decrease of the 
highest charge eigenvalue that can be traced to the vicinity of the Fermi surface, as shown ahead. 

In order to understand the spectrum in greater detail we consider the various eigenvalues in momentum space. 



In Fig. 10 we consider the system in the normal phase and in the superconducting phase for two equal density 
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Figure 12: (color online) Total fidelity as a function of momentum. Left panel: different temperatures and same finite A; right 
panel: different temperatures and A = 0. 



matrices p a = p^ and plot the three eigenvalues of the ^-fidelity operator as a function of momentum, for a 2d system. 
The fidelity is one for all momenta since Trp a = Trp5 = 1. The depression in the highest eigenvalue (left panels) 
marks clearly the Fermi surface. For momenta larger than the Fermi momentum there are no electrons (except for 
thermal excitations contained in the Fermi function). So the eigenvalue corresponding to empty states is 1. The other 
eigenvalues are close to zero outside the Fermi surface. Due to the particle-hole transformation of the Bogoliubov 
transformation, the empty site eigenvalue (corresponding to doubly occupied sites in terms of the electrons) is also 
close to 1 inside the Fermi surface. Accordingly, the other eigenvalues are also close to zero inside the Fermi surface. 
The noticeable features are therefore close to the Fermi surface. A similar structure is observed in the superconducting 
phase. There is a slight change close to the Fermi surface which is due to the opening of the superconducting gap. 
The amplitude of the doubly occupied and spin eigenvalues are smaller in this case since the decrease of the first 
eigenvalue is smaller along the Fermi surface. 



In Fig. 11 we show the fidelity operator spectrum for the case when one density matrix is in the normal phase and 
the other corresponds to a quantum state in the superconducting phase. The last panel shows the total fidelity. It 
is significantly decreased around the Fermi surface where the difference between the normal phase and the supercon- 
ducting phase is larger due to the pairing and opening of the gap. Note that the highest eigenvalue has a structure 
that strongly resembles the total fidelity. 

Finally, in Fig. [12] we compare the total fidelity as a function of momentum for different temperatures. The right 
panel corresponds to two quantum states in the normal phase but at different temperatures. As expected, around 
the Fermi surface the fidelity decreases however there is a sharp region where it approaches one. The width of the 
region around the Fermi surface is determined by the temperature through the Fermi function. The sharp maximum 
corresponds to the point where the two Fermi functions cross and so it pinpoints the location of the Fermi surface. 



CONCLUSIONS 



We have introduced and analysed the fidelity spectrum and the fidelity operator spectrum for partial states for 
different systems such as a magnetic impurity in a conventional superconductor, a XX spin- 1/2 chain in a transverse 
magnetic field and the thermal states of a finite temperature bulk superconductor. 

In the first problem we have found that only one charge eigenvalue and one spin eigenvalue have important changes 
as the quantum phase transition induced by the magnetic impurity occurs. The transition is associated with the 
capture of one electron by the impurity with a parallel spin. This feature is clearly seen when we consider two 
density matrices associated with the same lattice site (particularly the impurity site) and different but close by spin 
couplings between the impurity and the spin density of the electrons. Selecting two density matrices with the same 
spin coupling but different lattice sites leads to a signature of the phase transition that can be seen both from the 
eigenvalues associated with lattice sites far from each other and from the same lattice site. 

In the spin chain problem we have studied for the first time the block fidelity and the block fidelity susceptibility 
and found that the quantum phase transition that occurs between a XX phase and an Ising like phase is well signaled 
by the block fidelity. Both the entanglement spectrum and the fidelity spectrum do not show any significant features 
and we analysed the spectra calculating the moments of the distribution and the Renyi entropies. The Si Renyi 
entropy associated with the block fidelity shows a distinctive characteristic away from the critical point. 

Finally, in the finite temperature bulk superconductor we showed that in the superconductor there is a clear gap 
between the various /c-eigenvalues, as for the energy spectrum. In the case of two different density matrices we found 
that the effect of temperature is stronger than the difference in the order parameter distinguishing the normal from 
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the superconducting phase. Analysing the ^-fidelity operator spectrum it was clearly seen that the properties are 
determined by the structure around the Fermi energy, as expected. In the case of two density matrices for different 
temperatures in the normal phase the fidelity has a sharp maximum at the location of the Fermi surface, determined 
by the crossing of the Fermi functions. 

We have shown that the fidelity spectrum, which we have introduced, can give a more detailed description and 
characterization of the phase transitions of many-body quantum systems providing complementary information to 
other techniques. Therefore, we hope that this can be applied to non-trivial problems where the traditional Ginzburg- 
Landau theory with a local order parameter is not known. 
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